High energy neutrinos from radio-quiet AGNs 
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Most active galactic nuclei (AGN) lack prominent jets, and show modest radio emission and 
significant X-ray emission which arises mainly from the galactic core, very near from the central 
black hole. We use a quantitative scenario of such core-dominated radio-quiet AGN, which attributes 
a substantial fraction of the X-ray emission to the presence of abortive jets involving the collision 
of gas blobs in the core. Here we investigate the consequences of the acceleration of protons in 
the shocks from such collisions. We find that protons will be accelerated up to energies above the 
pion photoproduction threshold on both the X-rays and the UV photons from the accretion disk. 
The secondary charged pions decay, producing neutrinos. We predict significant fluxes of TeV- 
PeV neutrinos, and show that the AMANDA II detector is already constraining several important 
astrophysical parameters of these sources. Larger cubic kilometer detectors such as IceCube will be 
able to detect such neutrinos in less than one year of operation, or otherwise rule out this scenario. 

PACS numbers: 98.54.Cm, 95.85.Ry, 98.70.Sa 
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I. INTRODUCTION 

Around 90% of all active galactic nuclei (AGN), such as 
quasars and Seyfert galaxies, are radio-quiet, i.e. their ra- 
dio emission is relatively modest. This is due to the lack 
of prominent jets which provide most of the radio emis- 
sion in the remaining 10% of AGNs classified as radio- 
loud. Both types of AGN are also strong X-ray emittors, 
and the ultimate energy source is thought to derive from 
accretion of gas onto a central massive black hole (e.g. 
see 0). The high energy photon emission in radio-loud 
jet-dominated sources such as Blazars is largely due to 
non-thermal processes in the jet, whose dimensions ex- 
tend well beyond the outer edge of the galaxy. In radio- 
quiet AGNs, on the other hand, the short X-ray vari- 
ability timescales as well as spectral analysis of Fe X-ray 
lines indicate that the emission arises in the galactic core, 
very near from the central black hole. Both types of AGN 
have been considered as possible sources of high energy 
cosmic rays and neutrinos (see e.g. |2J and refs. therein). 
Here we concentrate on the more numerous radio-quiet 
AGN as a source of high energy neutrinos. 

The typical model of radio-quiet AGN X-ray emission 
is based on UV photons from the accretion disk suffering 
repeated inverse Compton scatterings in a hotter corona, 
leading to a power law X-ray spectrum (e.g. Q). This 
generic model does not allow for recent evidence indi- 
cating that, at some level, many radio-quiet AGN show 
weak radio emission, and in rare cases even show a weak 
incipient radio jet. A quantitative model which attempts 
to provide a bridge between the two classes of AGN con- 
siders the high energy emission in the core of radio-quiet 
AGNs due to the presence of an abortive jet Q|. These 
could arise from an outflow from the innner edge of the 
accretion disk, similar to that inferred in radio-loud ob- 
jects, but whose velocity is smaller than the escape ve- 



locity from the central black hole region. As in succesful 
jets, the outflow is expected to be intermittent, leading 
to shells or blobs ejected from the central region. In this 
case, due to their lower velocity, they will reach a maxi- 
mum distance from the hole before falling back, occasion- 
ally colliding with subsequently ejected, outgoing blobs. 
The X-ray emission in this model is attributed to Inverse 
Compton scattering of UV photons from the accretion 
disk by electrons accelerated in the sub-relativistic shocks 
resulting from blob-blob collisions. 

Here we consider the consequences of the acceleration 
of protons in the same shocks which accelerate the elec- 
trons. We find that protons will be accelerated up to 
energies above the energy threshold for pion photopro- 
duction on both the X-rays produced by the accelerated 
electrons, as well as on the UV photon field from the ac- 
cretion disk. The charged pions decay producing neutri- 
nos. A natural outcome of this model is then the produc- 
tion of a flux of high energy (TeV-PeV) neutrinos, which 
may be detectable with ~ km 3 volume detectors such as 
Icecube p|. Radio-quiet AGN constitute the main bulk 
of AGNs, and on average can be found closer to Earth. 
Also, since the neutrinos are produced in sub-relativistic 
shocks the emission is isotropic and the aborted jets need 
not point towards the Earth in order for the neutrinos to 
be detectable. We stress, however, that our conclusions 
are not restricted to this specific model of aborted jets, 
although it is used here as a framework. Qualitatively 
similar results might in principle be expected in disk- 
corona or other AGN models, if magnetic reconnection 
(e.g. M) or shocks associated with the accretion process 
(e.g. 0)13) l ea d to proton acceleration near the central 
black hole. 



In this paper we first show in Sec. HTlthat protons can 
be accelerated up to E™ ax « 10 9 GeV in the blob-blob 
collision shocks. We discuss proton interactions with X- 
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FIG. 1: Sketch of the geometry of a radio-quiet AGN core 
adopted in this work, showing the 10 8 Mq black hole (BH), 
accretion disk and two blobs moving in opposite directions 
which are about to collide. r a ~ 2.95 x 10 13 cm is a 
Schwarzschild radius for a 10 8 A/ n black hole. 



rays, UV photons and cold unaccelerated protons in Sec. 
IIIII In Sec. IIVI we discuss neutrino production by pho- 
tomcson and pp interactions, and we estimate the neu- 
trino flux from a typical radio-quiet quasar as well as 
the cumulative diffuse neutrino flux from all the radio- 
quiet quasars in the Universe. We also calculate in this 
section the muon and shower event rate expected in Gi- 
gaton Cherenkov detectors such as Icecube. In Seclvlwe 
summarize the paper and discuss the implications of our 
findings. 



II. PROTON ACCELERATION AND COOLING 
TIMES 

We adopt here a nominal set of parameters which serve 
as a specifc example for the purposes of normalization. 
We consider an AGN with a 10 8 solar mass black hole and 
with an accretion disk inner edge at a distance of about 
3r s from the black hole, where r s = 2GM/c 2 w 2.95 x 
10 13 cm is the Schwarzschild radius (Fig. We assume 
that blobs are ejected with an initial velocity j3o = 0.5, 
which corresponds to a Lorentz factor To ~ 1.15, i.e. they 
are non-relativistic and we assume that To ~ 1 from now 
on, and hence we do not make any distinction between 
the comoving (blob) frame and the observer's frame. The 
size of a blob is rbiob ~ ct r r s with a r ~ 1. The lifetime 
of the shock produced in a blob-blob collision is roughly 
equal to the blob-blob crossing time, t cr oss ~ a r r s /c/?o ~ 
2.0 x 10 3 s. 

The Eddington luminosity of a 10 8 M Q black hole AGN 
is L E dd = AnGMm p c/cFT « 1-26 x 10 46 L Eddfi erg s" 1 . 
This corresponds to an Eddington accretion rate M E dd = 
LEddhc 2 w 2.8M Q L Eddt8 7y^o 8 yr _1 , where 77 = 
0.08 770.08- Following Q| we assume that the accretion 
results in the intermittent ejection of blobs of matter. 
We also assume that a significant fraction aw, ~ 1 of 
the observed X-ray luminosity is due to the blob-blob 
collisions, whose number at any given time is JVbb- Abb 
can be estimated through a simple geometrical argument, 
from the number of blobs (A^biob) that are present in the 
region around the BH at that given time. Abiob is given 
roughly by the ratio of the time of flight of a blob i.e. 
the time it takes for a blob to go upwards from the BH 



and to fall back onto it, flight ~ 6r s //?oC, to the average 
time between the launch of consecutive blobs At. Taking 
At w (1 — 2) r s /c, a value that qualitatively accounts for 
the time variability of radio-quiet quasars we obtain 
-Wbiobs ~ 10. The average number of collisions (optical 
depth) of a blob during its flying time is given by the ratio 

1/3 

of the blob flying time tflight ~ /(A) c ) to the mean 
time between blob collisions, Vt /(NbiabVu^fioc), where 
the volume of a blob is Vbiob ~ (4/3)7rr 3 and the total 
volume inhabited by the blobs is Vt ~ (4/3)7r(3r s ) 3 (see 
Fig. ^) , giving an average number of collisions per blob 
w 1. The number of blob-blob collisions at any given 
time is then 1 < Nm < 5. We then assume that the 
X-ray luminosity from the steady-state number TVbb of 
blob-blob collisions is equal to a fraction abb ~ 1 of the 
total X-ray luminosity Lx- Assuming Lx to be a frac- 
tion ax, Ed < 1 of the Eddington luminosity, we combine 
these two factors into a single efficiency ax = Lx/L Ed d 
which we parametrize as cvx ~ 0.01 ao.01 corresponding 
to Lx ~ 10 44 erg s" 1 , which is the X-ray luminosity of 
a nominal radio-quiet AGN Q . All the calculations per- 
formed are parametrized by ax, which can vary in order 
to adjust Lx to the actual observed X-ray luminosity 
from a particular AGN. 

Using the usual assumption that the magnetic field B 
in the shock is due to turbulent fields generated behind 
the shocks arising from the blob-blob collisions, the co- 
moving (in the blob frame) magnetic field energy density 
is a fraction e B of the post-shock proton thermal energy: 



~ £ B PprcshockC 

07T 



(i) 



The density p pr eshock should be roughly equal to the den- 
sity of the accretion disk at a radius ~ 3r s , which is given 
by 0, 



Ppreshock ~ Pdisk (3r s ) 



M 



ATr(3r s ) 2 a v ^/GM/3r s 



w 1.5 x 10 14 a x ,o.oiL E dd,8 Vo.os § cm 

where M = axL E dd/v c2 an d we have taken the viscosity 
a v s» 0.1 0| • From p prcs hock the proton number density 
is, 

Pprcshock n 1A Q t — 1 —3 

n p = w 8.7 x 10 a x ,o.oiLEdd,s r] o os cm 

Tflp 

(3) 

where m p is the proton mass. The magnetic field is then 



B 



1/2 1/2 



/Q I 

° ne B ^preshock 

1.0 x 10 4 (ax,o.oi W 1/2£ fl a o.8 G ( 4 ) 



where e B = (1/3) es,o.3- 

The proton acceleration time is given by t° cc w A p r E /c 
where r E is the Larmor radius in the magnetic field in the 
blob, and A p is a constant > 1. Using Eq. Q we get, 
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A p ~f p ni p 

eB 
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ps 9.9 x 10" 9 ^ P (ax,o.oiW8)" 1/2 e B 'o.3 
The proton synchrotron cooling time is 



7 P <5) 



4-syn _ 

L p — 



6irm p c 
<JT h m 2 e ~j p B 2 

4.4 x 10 10 {ax.o.oiLEdd.^Y 



le Bfi.3 7 P 1 s (6) 



where CTh ~ 0.665 x 10 -24 cm 2 is the Thomson cross 
section. Equating the synchrotron cooling time and the 
acceleration time (see Eq. J5J) we get the maximum en- 
ergy of the accelerated protons, 



2.0xl0 9 A- 1/2 {a x ,o.oiL E dd,8 



,-1/4,-1/4 



',0.3 



GeV 
(7) 

Protons can also suffer Inverse Compton (IC) losses 
from interactions with X-ray and UV photons from the 
disk. To estimate the X-ray photon number density, we 
take into account that the blobs are optically thin |4| so 
X-ray photons escape freely, and their density is nx — 
axL Edd /4ir(3r s ) 2 cEx, or 



nx 



CtxLEddt cross 

(4/3)irrlEx 



ps 5.3 x 10 13 axM.oiL E dd,8 £ x ,keV cm 3 ( 8 ) 

where Ex = 1 £x,kcV keV is the typical energy of the 
emitted X-rays. 

To estimate the density of UV photons from the ac- 
cretion disk, we consider thermal emission in the in- 
ner disk between pa 3r s and 7.5r s . This portion of 
the disk emits an UV luminosity a\jyLEdd- Observa- 
tionally the X-ray luminosity is of the order of 10-50% 
the luminosity in the optical-UV We then take 

auv ~ 5ox ~ 0.05 auv,o.05- The equilibrium temper- 
ature (energy) of the UV photons can be approximately 
estimated assuming the emission from the disk follows a 
black-body spectrum, yielding: 

1/4 



UV 



/ auvL 



Edd 



\ a Adisk 

« 8.3 (auvAOS^B^s) 174 eV (9) 

where ks and a are Boltzmann's and Stefan-Boltzmann's 
constant respectively, and Adisk is the area of the inner 
disk. The number density of UV photons is then, 

OlwLEdd 



n uv 



^diskCi^UV 

1.2 x 10 16 (a uv ,o.5i™) 3/4 



(10) 



and we get that nuv ~ 200 nx- 

The inverse Compton scattering of protons on X-rays 
and UV photons is in the Thomson limit for energies 
below Ejp, where 
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Ex,uv 

f 8.8 x 10 5 E 
1 1.1 x 10 



X.keV 
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FIG. 2: Proton cooling times as a function of proton energy 
for different energy loss processes, namely, synchrotron and 
inverse Compton scattering on both X-rays and UV photons. 
Also shown are the proton acceleration time, and the blob- 
blob crossing time. 



The IC cooling times on X-ray and UV photons are then, 
in the Thomson (Th) limit 



t 



IC,Th 



m„ 



<JTh(rn e /m p ) 2 c~f p Ex,Tjv n x ,vv 
f 3.0 x 10 12 {a x ,o.oiL E dd,&) 
~ \ 1.6 x 10 12 {auYflmLEddfi)' 1 !^ 1 s 

and in the Klein-Nishina (KN) (E p » E l p c ) limit 

7 P (£; 7 /m p ) 



(X) 
(UV) 
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<JTh(m e /m p ) 2 c n x ,uv 
3-4 (ax,o.oi^£;drf,8)~ 1 £'x,kcv7p 
1.2 x 10- 4 (a vv , .o 5 L E dd,8)- 1/2 
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In Fig. |21 we show the IC cooling times on X-rays and 
UV photons along with the synchrotron cooling and the 
proton acceleration times as a function of proton energy. 
Also shown is the dynamic (blob-blob crossing) time. It 
can be seen that IC losses and adiabatic losses (given by 
the crossing time) do not impose any constraints on the 
maximum proton energy, which is dominantly restricted 
by the synchrotron energy losses, and is given in Eq. J7J. 



III. PROTON INTERACTIONS 

High energy protons accelerated in the blob-blob col- 
lision shock can interact with X-rays produced in the 
shocks, UV photons in the disk or cold (unaccelerated) 
protons in the blobs. The proton threshold energy for 
photomeson production on X-rays and UV photons at 
the A-resonance is given by E p E 7 ps 0.3 GeV 2 which 
yields, 

Ef ith ps 3.0 x 10 5 E-^ v GeV (11) 
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for photomeson production on X-rays, and 

Epjh « 3-6 x 10 7 (a lJv .o.o5L E d d ,8r 1/4 GeV (12) 

for interactions on UV photons. 

The threshold for inelastic pp — > pnn + interactions in 
which neutrinos can be produced is much smaller and is 
given by, 



E p - 



(m p + m n + TOtt) 2 — 2m 2 
2m„ 



1.23 GeV (13) 



The optical depth for photomeson interactions on X- 
rays, from Eq. (JSJ, is 

T p-t ~ r s(o"p 7 -> A)n x « 7.8 x 10 _1 ax.o.oi^Edd.s^x keV 

(14) 

and the corresponding optical depth on UV photons, 
from Eq. ifTUI) . is, 

r p u 7 v « r s (cr p7 -> A)n uv » 1.8 X 10 2 (a VVfi mL E dd,8) 3/4 

(15) 

where <t P7 ^a ~ 5 x 10 28 cm 2 is the photomeson pro- 
duction cross section at the A-resonance. 

Protons can also interact with cold (unaccelerated) 
protons in the blob. We assume that a fraction £ p < 1 
of the protons are Fermi accelerated in the shocks pro- 
duced in the blob-blob collisions. The pp collisions should 
dominate below E* th> the energy threshold for pho- 
tomeson production on X-rays. In this energy range the 
mean total cross-section for pp interactions is ((J pp } ~ 
6 x 10 ~ 26 cm 2 , roughly two orders of magnitude larger 
than the photomeson production cross section and the 
density of proton targets is typically 10 3 times smaller 
than the X-ray density, resulting from Eq. J2J in a 
smaller optical depth for the Lx 
here, 



10 44 erg s 1 adopted 



Vpp) ~ ( f7 'pp)^s(l £ P )^p 

« 1.5 x 10~ 2 (1 - £ p ) ax, o. oi L Edd , 8 (16) 

This is the optical thickness for accelerated protons en- 
countering blob thermal nucleons. In addition, acceler- 
ated protons emitted in some directions will also col- 
lide with inner disk thermal nucleons. If the disk in- 
ner edge is at 3r s and has a thickness 3r s , a fraction 
27r(3r s ) 2 /(47r(3r s ) 2 ) ~ (1/2) will collide with nucleons in 
the inside wall of the accretion disk, for which r pp is larger 
due to the larger path length available. The density at 
inner edge is comparable to the blob density, and varies 
with radius, while the effective path length depends on 
the incidence angle. We adopt as an estimate for the 
average proton incident on the disk a value of r pp ~ 10 
times larger than in a blob. 

Inverse Compton scattering of protons on X-rays or 
UV photons does not affect the pry scattering rate, since 
the optical depth for IC scattering tic < 1CU 1 r p7 for 
all target photons and in both the Thomson and Klein- 
Nishina regimes. 



The ratio of IC to pp optical depth in the Thomson 
regime (E p < E l p c ) is 



tic _ QTh "x,uv (me/rrip) 2 
(rpp) (cr pp ) (1 - £ p )n p 
/ 2.0xl0- a (l-&)- 1 *£i 
j 4.6 (Wp)- 1 *; 



(17) 



X,kcV 



(X) 



1 —1 3/4 j— ±/i (TT\r\ 
a X,0.01 a UV,0.05^Edd,8 \ V V I 

and hence pp interactions which are dominant below the 
energy threshold for photomeson production on X-rays, 
are suppressed by IC scattering on UV photons for the 
Lx ~ 10 44 erg s _1 and Luv ~ 5 x 10 44 erg s _1 of the 
nominal radio-quiet quasar. In the Klein-Nishina regime 
T ic/( T PP ) is lower than above, however in this regime pry 
interactions dominate as we will see below. 
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A. p7 interactions on a power-law X-ray spectrum 

High energy neutrinos are produced in pry interac- 
tions with both X-rays and UV photons above the cor- 
responding energy thresholds E* th s» 3.0 x 10 5 GeV 
and EzY th m 3.6 x 10 7 GeV respectively. However, the 
X-ray spectra of AGN, including the radio-quiet, are typ- 
ically of power-law form, and hence protons at energies 
below th can also produce delta resonances by in- 
teracting with photons from the high energy tail of the 
X-ray spectrum. The same is true for pry interactions 
with UV photons in which low energy protons can pho- 
toproducc on photons from the high energy tail of the 
blackbody photon distribution. As a first approximation, 
we will neglect the tail of the blackbody distribution, and 
assume that the UV spectrum is monoenergetic. How- 
ever, we will take into account that the X-ray spectrum 
follows a power-law given by dNj/dE^ cx E~ s . As a 
consequence the optical depth to p- X-rays depends on 
the spectral index S of the photon spectrum, as well as 
on proton energy. This is due to the fact that as E p 
increases, protons interact with an increasing number of 
lower energy photons. For a steep X-ray spectrum, the 
bulk of the X-ray photons are concentrated in the low 
energy region of the power-law distribution, and this in- 
creases the optical depth to p- X-ray interactions for high 
energy protons. For a photon number density spectrum 
oc E- 5 , the energy dependent optical depth is given by, 



r* (£ 7 ) ~ r s (a pX ^A)E^ x E^ s 



(18) 



Using the threshold relation £ p £ 7 w 0.3 GeV 2 Eq. JTSJ 
turns into, 



,(£?„) c r p x 7 (£ 7 = lkeV)(-A 



p.th t 



7.8 x 10 1 (a X:0m L Edd ^) 



E n 



E x 



6-1 

I (19) 



5 



where t p \{E 1 



1 keV) 



X.kcV 
PI 



„.,,—. - • ,, is the optical depth 

given in Eq. The factor (E p /E p K th }^- 1 '> accounts 

for the increasing optical depth as E p increases (note that 
this is only for 8 > 1, which is observationally the case). 

From the discussion above, it follows that (for 8 > 1) 
there is a proton energy E p above which photoproduction 
on X-rays dominates over photoproduction on UV pho- 
tons. The energy E p is given by the condition T^{E p ) = 
. Solving for E p we have 

1/(5-1) 



PI 

X.kcV 

T p ~) 



^p,th 



[22Q 6 a" 1 a 3/4 r -1 / 4 



(20) 
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The predominant target for pj interactions is determined 
by the numerical value of E p and whether it is larger or 
smaller than E^ h , the energy above which p's photopro- 
duce on UV photons, as well as by its relative value with 
respect to E r p nax , the maximum energy of the accelerated 
protons. E p depends on the X-ray spectral index 8. It 
follows that: 

• If 8 > 8 C w 2.13 then E° < E^J h , implying 
that proton interactions with X-rays photons al- 
ways dominate over interactions with UV photons. 



If 1 < 8 < 8, 



1.61 then El 



> E p nax and inter- 



actions with X-rays dominate up to E^ th , while 
photoproduction on UV photons dominates above 
E uv 

Pl,th 



• If 8i < 8 < 8 C there is a small proton energy range 
above E^Y th but below E p nax in which photopro- 
duction on UV photons dominates, whereas inter- 
actions with X-rays are more important outside it. 



The numerical value of 8 C is obtained from the condition 
Ep(fic) = EpY h whereas 8i is obtained from the equation 
E p (Si) — E™ ax . They depend on the particular values 
of Lx and L\jy. The situation is illustrated in Fig. [31 
where E p is plotted as a function of 8, along with E™^ 
and E p ^Y th . The dark grey shaded area represents the re- 
gion in the space parameter of {E p , 8) where interactions 
on X-rays are dominant, while photoproduction on UV 
photons is more important in the light gray shaded area. 
Interactions above E r p nax (unshaded region) are not pos- 
sible, since protons do not get accelerated to E p > E™ ax . 

For the pp interactions, similarly to the discussion 
above, there is a proton energy EP p below which pp in- 
teractions are more important than interactions with X- 
rays, given by the condition r p , i(-E£ p ) = (t pp ). Solving 
for EP p we get: 



'pp/ 



X.kcV 



' in 



1/(5-1) 



TP* 
P,th 



[9.8 x l(T 3 ]W-i) E x th 



(21) 



FIG. 3: Several energy scales relevant for proton interactions 
for the set of parameters of the nominal radio-quiet quasar: 
Mbh = 1O 8 M , ax, o.oi = 1, auv.o.05 = 1 



E™ ax is the 



maximum energy for proton acceleration. E^ th is the pro- 
ton energy threshold for photomeson production on UV pho- 
tons. £™ n is the lowest energy of a proton that can interact 
through the A— resonance with an X-ray assuming the maxi- 
mum X-ray energy is E™ ax = 100 keV. The upper curve is Ep, 
the energy above which the optical depth for pry interactions 
on X-rays is larger than the corresponding optical depth on 
UV photons. The lower curve is Ep P the energy below which 
the pp optical depth is larger than the p7 optical depth on X- 
rays. p7 interactions on X-rays are dominant in the darkest 
grey shaded parameter space region (E p ,8). p7 interactions 
on UV photons dominate in the lightest gray shaded area, 
whereas pp interactions are more important in the medium 
gray shaded region at the bottom of the plot. 



There is also another important energy scale, E™ ax 
which is the photon energy above which the source runs 
out of energetic enough X-rays to interact with protons of 
energy below E p nm given by the photoproduction thresh- 
old condition: E p nin E^ ax » 0.3 GeV 2 . Observationally, 
most radio-quiet AGN spectra show an X-ray cut-off at 
about 100-150 keV. We will adopt E™ ax = 100 keV in 
this work. 

As in the case of j>7 interactions, E pp depends on the 
index of the X-ray spectrum. For a steep enough spec- 
trum there are not enough high energy photons to in- 
teract with low energy protons and then pp interactions 
dominate below a proton energy larger than E p nm . This 
occurs for 8 > S p ~ 2.00 given by Epp(8 p ) = E™'" 1 , using 
£ p = 0.5 and the standard set of luminosities of the nom- 
inal radio-quiet AGN. This is also shown in Fig. [3]where 
Epp is plotted along with E r p nm assuming E™ ax = 100 
keV. Thus, pp interactions dominate in the medium gray 
region at the bottom of Fig. [3] if the p flux is not sup- 
pressed by IC scattering, which depends on the luminos- 
ity of the source. 
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IV. NEUTRINO PRODUCTION AND 
DETECTION 

High energy neutrinos are produced through secondary 
charged it decay photoproduced in pry — > A and 777 — > 
A interactions. Each of the 4 leptons produced in the 
charged pion decay chain — > /ii^, — > eu e v^9p carry 
about 1 /4 of the initial pion energy. 

Neutrinos can also be produced through pion and 
muon decay produced in pp interactions in the shocks. 
In fact pp interactions are dominant at low proton ener- 
gies (medium-grey shaded region at the bottom of FigJSJ . 
The typical neutrino energy is 1 /4 of the initial pion en- 
ergy, and this is determined by the pion multiplicity in 
the pp collisions. 

Due to the high magnetic field and dense photon en- 
vironment where the 7r + and [i + are produced, they 
might lose their energies through synchrotron and inverse 
Compton scattering before decaying to neutrinos. 

The synchrotron cooling time for pions and muons 
is calculated using Eq. JBJ. The charged pion and 
muon decay times are r^ ec ps 2.6 x 7^ s, and 

T dec w 2.2 x 10~ 6 7 M s. Equating the synchrotron cooling 
time and the decay time, we get the synchrotron break 
energies above which neutrino production is suppressed 
by a factor oc E~ 2 due to energy losses: 



Ei 



1.0 x 10 7 (a x ,o.oi^Ed<i,8)" 1/2 e B 1 o /2 3 GcV ( 22 ) 



for pions, and, 

Ef « 5.6 x 10 5 (ax, o.oi£ 



EddA 



~ 1/2 ^ll GeV (23) 



for muons. 

Inverse Compton scattering on both X-rays and UV 
photons does not contribute significantly to the pion and 
muon energy loss. For all photon targets and in both the 
Thomson and Klein-Nishina regimes, the inverse Comp- 
ton cooling time is at least ps 10 times larger than the 
decay times (and much larger than this ps 10 2 — 10 7 in 
most cases) and hence we ignore pion and muon IC scat- 
tering energy losses in the calculation of the neutrino 
flux. 



A. Neutrino flux from the nominal radio-quiet 
quasar 

Protons are accelerated by the Fermi mechanism in 
the blob-blob shocks up to energies E™ ax ss 10 9 GeV, 
following an energy spectrum oc E~ 2 . The energy per 
unit time given to protons producing neutrinos is, 



1 TFT~7T dE p = £ p CtyJjEdd 

J E p dEpdt 

p,th 

Taking d 2 N p /dE p dt = K p /E 2 we get for K p , 



(24) 



w 3 2 x 1Q 45 ^X, .0lL E dd,S ! ( } 

HE^/E^) K > 



The proton flucncc per unit time is, 



F„ 



1 



Et 



d 2 N 



4irD 2 p dE v dt 



(26) 



1.6 x 10~ 5 



£p &X,0.0lLEdd,t 



GeV 



D 2 L1 \n(E™ ax I 'E v v th ) cm 2 s 
for a source at a luminosity distance Dl = 1 Dl,i Mpc. 

1. Neutrino flux from P7 interactions 
The neutrino flux for each neutrino flavor is given by, 



$pt = -f F 



(27) 



e; 2 

E- 2 (E?/E v f 



Ef <E V <E™ 



sb 



where Ef, b ps (1/4) E^ b for muon neutrinos produced in 
~- ' 1 11 r ' b for electron and muon 



ir decays, and E™ « (1/4) 



neutrinos produced in the decay of /j's. E^> and E^ 
are given in Eqs. i(2*2|) and (|2*3^l respectively, E v vth = 
0.05 E p p th and E™ ax = 0.05 E™ ax . 

f n is the fraction of proton energy given to pions. As 
discussed in subsection IIII Al f n depends on proton en- 
ergy (and hence on neutrino energy) and on S the spectral 
index of the X-ray energy distribution, which determine 
the dominant photon target (X-rays or UV photons) . For 
interactions on X-rays f„ is given by, 



J* = mm 



1, (p^p — >7T") 



_X : kcV 
PI 



Ev 
E u,th J 



5-1' 



(28) 



where (x 



p^TTl 



0.2 is the average inelasticity in a pry 
interaction, r^, ,kcV is given in Eq. (fT4"jl and E* th = 
0.05 E* tfl given in Eq. (jTT|) . For interactions on UV 
photons, 



/" V = min(l, (xp^TpV) 



(29) 



where t^Y is given in Eq. (|T51) . 

The dependence of /„■ on 5 and E v can be summarized 
as follows, 

• J£8>S C , then: 

/* = /*; K, th < E »<K ax (so) 

• If 1 < S < 61, then: 

(/*; E v <EYJ h 
U = / uv ; E% h <E v <El (31) 

where E^J h ps 0.05 E^ th given in Eq. Q2J}, and 
El ps 0.05 E c p given in Eq. JUJ. 

• If Si < 5 < S c , then: 
/*; E V <E™ 



(32) 
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2. Neutrino flux from pp interactions 



The neutrino flux from pp interactions is 



f PP ,x M V (E P ) 



dE„ 



E 2 p 



\ J f PP ,2 M v {E p ) |j dE p 



(33) 
(34) 



for E. 



< E p,th> where fppA 



min(l, (t pp )), and / PP)2 



min(l, 10{t pp )) accounting for the larger optical depth in 
the inner wall of the accretion disk where ~ 1/2 of the 
accelerated protons interact (see Section ITTljl . M v (E p ) is 
the neutrino multiplicity through pion decay in the pp 
interactions which can be parameterized as [lJl , 



MJE P ) = - ( -^-) 



1 /lO^GeVN 

2 \ E p J 







1 TXlm 



4GeV 



7cm < 



E u 



< 



1 EL 



GeV ~ 4 GeV 



(35) 



for each neutrino flavor. 

Figure 0] shows the neutrino flux produced in pp and 
P7 interactions, for the nominal radio-quiet quasar at a 
distance D L = 20 Mpc. The quasar X-ray spectrum is 
fairly universal, typically being E~ s with 8 = 1.7 [l4| . 
For this spectral index and for the nominal radio-quiet 
quasar, interactions with unaccelerated protons, X-rays 
and UV photons contribute to the neutrino flux (see Fig. 

Below the proton threshold for interaction on X-rays, 
pp interactions dominate neutrino production, however 
pp collisions are largely suppressed due to the more dom- 
inant inverse Compton scattering on UV photons (see 
Eq. I18|) . This behavior continues until the protons have 
enough energy to interact with the high energy tail of 
the X-ray photon spectrum. This effect accounts for the 
rapid increase in the neutrino flux around « 200 GeV 
seen in Fig. 0] After that the flux is slightly flatter than 
the expected E~ 2 power law due to the -E^ -1 = E^ 1 
factor in the optical depth. As energy increases f n 
reaches 1 and the neutrino flux goes as E~ 2 , before syn- 
chrotron losses by pions become important making the 
neutrino flux to fall as E~ A . 

The normalization as well as the shape of the neu- 
trino flux changes with luminosity. This is illustrated in 
Fig. 0] for a radio-quiet quasar of Lx ~ 10 45 ergs -1 
and Mbh = 10 8 M . For larger Lx the optical depth 
to pX— rays interactions reaches 1 at a smaller proton 
(and hence neutrino) energy due to the denser photon 
target, and as a consequence the E~ 2 behavior of the 
spectrum starts at a lower energy than in the nominal 
AGN. Besides, since for fixed black hole mass, the pion 
synchrotron break energies scales as L^ 1 ^ 2 (see Eq. (I22f) 'l 
the flux starts to fall as E~ 4 at a lower energy than for 
the lower luminosity AGN. 



s 

> 



L x =10 ergs 
L x =10 ergs 




l(f 1(T 10' 10 J 10" 10' 10° 

[GeV] 

FIG. 4: Solid line: Neutrino flux reaching Earth for the 
nominal radio-quiet quasar: Mbh = 1O 8 M0, ax = 0.01 
(Lx ~ 10 44 erg s _1 ) and auv = 0.05 at a luminosity distance 
Dl = 20 Mpc, (l+z) ~ 1. Dashed line: Neutrino flux for a set 
of parameters M B a = 10 8 M Q , a x = 0.1 (L x ~ 10 45 erg s" 1 ) 
and auv = 0.5 at Dl = 20 Mpc. Only from pion decay 
produced in pp and py interactions are shown in the figure. 



B. Diffuse neutrino flux 

The cumulative neutrino flux from all radio-quiet 
AGNs in the Universe is obtained by convolution of the 
observed point source neutrino flux with the luminosity 
function taking into account its cosmological evolution, 



1 

4?r 



ciiir 



ob(^,ob) — 
L max 

dL x 

£mi„ 



dz^f(z)^f^ ob (36) 
dLx dz 



where dn Q /dLx describes the present day X-ray lu- 
minosity function of the sources and f(z) is its 
cosmological evolution. We have used the broken 
power-law luminosity and evolution functions of refer- 
ence [3 (model I). In a Friedmann-Robertson- Walker 
universe the comoving volume element is dV/dz — 
ATiD\c\dt/ 'dz\j '(1 + z) and the derivative of the cos- 
mic time t with respect to redshift z is (dt/dz)^ 1 = 
-H (l + z)y/{\ + VL m z)(l + z) 2 - fl A (2z + z 2 ) where we 
have used a standard ACDM cosmology with f2 TO = 0.3 
and £!a = 0.7. The neutrino energy in the observer's 
frame and the source frame are related by E v = £'„ i ob(l + 
z), and the luminosity at the source is (1 + z) 2 the lu- 
minosity observed today. In the integration we take into 
account that the shape of the observed individual neu- 
trino spectrum depends on the luminosity of the source 
as illustrated in Fig. 0] 

We have performed the calculation of the diffuse flux 
in two different ways. First we assume that in all the 
individual sources the same fraction of the Eddington 
luminosity (ax = 0.01) is converted to X-rays, adjusting 
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the mass of the black hole in order to get the necessary 
X-ray luminosity (variable Mbh case) . As an alternative 
calculation, we assume the black hole has the same mass 
in all sources (Mbh = 10 8 M©) and we vary the efficiency 
of conversion of Eddington luminosity into X-rays, i.e. we 
vary ax, adjusting its value so that we get the required 
Lx (fixed Mbh case). 

In Fig. [5] we plot the diffuse + flux as obtained 
in both calculations. The contribution to the total neu- 
trino flux from radio-quiet AGNs in different luminosity 
bins is shown. The less luminous AGNs, although more 
abundant, do not contribute much to the total because 
the neutrino flux scales with luminosity. The most lumi- 
nous sources, although powerful neutrino emitters, are 
less abundant and contribute very little to the total flux. 

In the case of fixed black hole mass (bottom panel) the 
dimensions of the accretion disk and black hole region 
are also fixed, and as luminosity increases the density of 
protons and X-rays increases linearly with Lx- Proton 
interactions with X-rays start to become important at 
an increasingly smaller energy as Lx increases due to the 
denser X-ray target. At the same time, the synchrotron 
break occurs at an increasingly lower energy because the 
magnetic field increases as L^ 2 (see Eq. Q). The two 
effects are visible in Fig. [S] 

In the case of variable Mbh (top panel) , an increase in 
luminosity for a fixed accretion efficiency is accompanied 
by an increase of the black hole mass, i.e. the accretion 
disk and black hole region becomes larger, and the proton 
and photon density actually decrease as Mgjj. This has 
the effect of decreasing the magnetic field intensity in the 
source and hence the synchrotron break occurs at larger 
energy as Lx increases. However the optical depth to 
P7 interactions stays the same as Mbh increases and as 
a consequence the proton energy at which it reaches 1, 
i.e. the energy at which the neutrino flux starts falling 
as E~ 2 , stays the same. The two effects are again visible 
in the top panel of Fig. |3J 

Using Eq. 119|) one sees that for low luminosity AGNs, 
Lx < fO 41 ergs" 1 , < 1 for all proton energies up 
to E™ ax , and protons can escape from the sources con- 
tributing to the observed cosmic ray spectrum. As a 
consequence, neutrinos produced in the blob-blob shocks 
in low luminosity AGNs would in principle be affected 
by the Waxman-Bahcall bound: E 2 dN v /dE u < (I -4) x 
1CP 8 GcV cm~ 2 s^ 1 However, the contribution of 

sources with Lx < 10 41 erg s _1 to the total neutrino flux 
is in any case below the Waxman-Bahcall bound, as can 
be seen in Fig. while for larger luminosities, the neu- 
trinos are produced in regions where the optical depth to 
P7 collision is typically > 1. Thus the Waxman-Bahcall 
bound does not apply to the sources considered here. 

Also shown in Fig. 0are the latest AMANDA II limits 
on + for an E~ 2 diffuse flux which apply in 

different energy ranges, as well as the expected IceCube 
sensitivity (90% C.L.) after 1 year of operation for an 
E~ 2 flux The energy range in which the IceCube 
sensitivity is shown, is for illustrative purposes only. 



It is important to remark that the nominal muon event 
rates (see below) , as well as the individual source and dif- 
fuse neutrino fluxes of figures 0] and are proportional 
to (the fraction of the X-ray flux in radio-quiet AGN 
which is due to a core component) and to £ p (the frac- 
tion of protons which is accelerated in core shocks of any 
origin, attributed here more specifically to blob-blob col- 
lisions). The values adopted here as an example, abb ~ 1 
and £ p ~ 0.5, are reasonable but not well constrained so 
far. Additional uncertainties for the diffuse neutrino flux 
are the approximate luminosity function, and the fraction 
of radio-quiet AGNs (which may be thought of as an av- 
erage abb) for which a core X-ray component requires 
invoking shocks, blob collisions, aborted jets or other 
mechanisms besides an accretion disk. The AMANDA 
II neutrino detector is currently sensitive to the diffuse 
neutrino fluxes predicted here with the above parame- 
ters (see Fig. |5J, and hence gives important informa- 
tion about the relevant astrophysical parameters of the 
aborted jets scenario, namely, that the product abb x £p 
may be a factor < 1/2 times the nominal values adopted 
here. The larger km 3 IceCube will be able to detect the 
diffuse fluxes in less than 1 year of operation, or otherwise 
rule out the generic model explored in this paper. 



C. Event rate in a Gigaton telescope 

We estimate here the neutrino event rate from a sin- 
gle nearby radio-quiet AGN, taking as an example a dis- 
tance Dl « 20 Mpc, comparable to that of the Virgo 
cluster of galaxies. We assume the nominal radio-quiet 
AGN with parameters: Mbh = 1O 8 M , auv,o.05 = 1, 
and ax, o.oi = 1, i- e - Lx ~ 10 44 erg s _1 . Consider- 
ing neutrinos from p7 interactions, the energy in neu- 
trinos reaching the Earth per unit area and time from 
a source at this distance, from Eqs. (|2*7|) and is 
F u w 10 -8 GeV cm -2 s~ 4 , where we have used £ p = 0.5, 
/„ = 1, Ep lax = 2 x 10 9 GeV and E p p th = 1.23 GeV. 
The number of neutrinos reaching the Earth assuming 
a typical neutrino energy E v = 100 TeV = 10 5 GeV is 
« F u /E u w 10~ 13 cm~ 2 s~ 4 . At this energy the prob- 
ability that a neutrino converts into a muon within the 
range of the detector is w 10~ 4 . The number of de- 
tectable muons is rs 3 per km 2 per year. This number 
scales roughly with the luminosity of the source and the 
inverse of the distance squared. The estimate only takes 
into account muon events from i/^'s produced directly in 
7r decay. We have also made a more detailed estimate 
that convolutes the neutrino flux with the probability of 
neutrino detection through both muon and shower obser- 
vation, taking into account the absorption of the neutrino 
flux while passing through the Earth (see for instance 
Appendix C in |19j). and assuming that after propaga- 
tion from the source neutrinos have oscillated so that the 
flavor ratio is v e : : v T = 1 : 1 : 1. We obtain a 
muon event rate m 9 km~ 2 yr _1 and a shower rate 
w 2 km -2 yr~ 4 for an AGN located at a zenith an- 
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FIG. 5: Radio-quiet AGN diffuse + flux reaching Earth 
(thick solid line). The fact that half of the muon neutrinos 
disappear due to oscillations is accounted for. The contribu- 
tion to the total neutrino flux (thick lines) from radio-quiet 
AGNs in different X-ray luminosity bins (in erg s _1 ) is shown 
in thin lines. The top panel shows the results of a calculation 
where the black hole mass is varied, while the bottom panel 
shows the case where the accretion rate is varied (see text for 
details). Also shown are the latest AMANDA II limits on an 
E~ 2 diffuse u^ + u^ flux which apply in different energy ranges 
(thick dashed-dotted lines) |l5lll6f. and the expected IceCube 
sensitivity after 1 year of operation (for an E~ 2 diffuse v^+v^ 
flux at 90% C.L.) (thick dashed line). 



gle 90 degrees (with respect to the vertical direction) in a 
detector such as Icecube at an average depth of ~ 1.8 km 
water equivalent. The muon and shower energy threshold 
we used is 500 GeV. 



For the diffuse fluxes plotted in Fig. [5] we predict a 
muon rate integrated over the 27r northern sky (where 
the background due to atmospheric muons is negligible) 
of the order of a few hundred events and a shower rate 
about 3 times smaller, with roughly the same rate for the 
fixed Mbh case. IceCube can also identify events com- 
ing from the southern hemisphere despite the large muon 
background by selecting events with high enough energy 
(jj |2(j. The number of muon events expected in this 
model from the 2tt southern sky for an energy threshold 
of 10 6 GeV is AT w 17 km" 2 yr" 1 for the variable M BH 



case and N„ 



km 2 yr 1 for the fixed Mbh case. 



V. SUMMARY AND CONCLUSIONS 

We have investigated the high energy neutrino emis- 
sion from radio-quiet AGNs, which is the most numerous 
class of AGNs. The neutrino emission in these AGNs, 
which lack prominent jets, is attributed to pj and pp 
collisions of relativistic protons accelerated in shocks oc- 
curring near the central black hole, which are also re- 
sponsible for a core component of the observed X-ray lu- 
minosity. The calculations were performed for a specific 
model of the AGN core in which an abortive jet leads 
to collisions of blobs inside the inner edge of the accre- 
tion disk. The results, however, should be generic to any 
model involving particle acceleration in the inner region 
of the accretion flow, e.g. 0,0. This model allows us to 
predict spectral details whose features depend on physi- 
cal parameters related to the black hole mass, accretion 
rate and X-ray/UV luminosity, as well as shock physics 
parameters such as proton injection fraction, shock mag- 
netic field strengths, etc. In particular, the model pre- 
dicts muon and shower event rates from individual AGN, 
as a function of their observed X-ray/UV luminosity. The 
individual source event rates, as well as the correspond- 
ing diffuse fluxes are, for nominal parameters, already 
constrainable with AMANDA II, and should be strongly 
constrainable with the IceCube detector. 
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